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Abstract 

A  formalism  is  presented  which  demonstrates  that  the  mean  Green’s  function  for  the 
acoustic  field  between  two  rough  interfaces  can  be  expressed  as  a  Green’s  function  associated 
with  two  flat  interfaces  with  effective  reflection  coefficients.  This  result  incorporates  all  orders 
of  the  fluctuations  in  the  half-space  scattering  amplitudes  associated  with  each  interface 
considered  separately.  From  the  mean  Green’s  function  modal  attenuations  can  be  found. 
To  lowest  order  in  the  surface  height  fluctuations  it  is  shown  that  it  is  not  sufficient  to 
use  mean  half-space  scattering  amplitudes  as  effective  reflection  coefficients.  The  formalism 
is  designed  to  provide  approximations  for  the  Green’s  function  in  layered  media  which  are 
based  on  previously  developed  approximations  for  half-space  scattering  amplitudes. 


■  IIS  s  I  I  ■■psw— n— a— i— ^ 

This  documsnt  bos  b««a  appiofwd 
for  public  taiaosa  emd  sola;  itt 
disTOUtioo  is  ualimiiad 


#4  t  19  092 


94-22639  ^ 


one  QUALITY  DiSFELhlkL  i 


I  Introduction 


A  number  of  good  approximations  have  recently  been  developed  for  plane- wave,  half-space  rough 
surface  scattering  amplitudes.  The  aim  of  this  paper  is  to  show  how  these  scattering  amplitudes 
can  be  used  to  describe  sound  propagation  between  two  rough  interface.  This  work  is  motivated 
by  a  certain  laziness:  it  shouldn’t  be  necessary  to  repeat  the  discussions  used  to  describe  the 
scattering  from  one  interface  when  treating  two  interfaces.  One  should  only  be  required  to  treat 
the  interactions  between  the  interfaces. 

In  an  earlier  paper  [1]  the  method  of  smoothing  was  used  to  treat  this  problem  to  lowest 
order  in  scattering  amplitude  fluctuations.  It  was  indeed  shown  there  that  there  are  efltective 
reflection  coefficients  for  the  mean  Green’s  functions  which  are  distinct  from  the  mean  half¬ 
space  scattering  amplitudes.  Here  the  work  of  Ref.[l]  is  extended  to  all  orders  in  scattering 
amplitude  fluctuations.  It  is  important  here  to  distinguish  scattering  amplitude  fluctuations 
from  fluctuations  in  surface  heights.  In  general,  scattering  amplitudes  are  non-linear  functionals 
of  the  surface  height.  The  idea  here  and  in  Ref.[l]  is  to  develop  a  formalism  which  uses  non- 
perturbative  approximations  of  the  scattering  amplitudes. 

It  .  will  be  shown  here  that  the  existenc*  of  effective  reflection  coefficients  for  propagation 
between  two  rough  interfaces  is  not  an  artifact  of  the  smoothing  approximation.  In  fact,  the 
effective  reflection  coefficients  can  be  expressed  in  terms  of  self-enerpes  associated  with  a  cer¬ 
tain  random  operator.  In  constructing  the  effective  reflection  coefficients,  a  systematic  way  of 
including  higher  order  moments  of  the  fluctuations  of  the  half-space  scattering  amplitudes  will 
be  developed.  In  another  extension  of  the  work  in  Ref.[l],  it  will  be  shown  how  correlations 
between  the  scattering  on  the  upper  and  lower  interfaces  affect  the  mean  fteld.  If  there  is  corre¬ 
lation,  the  mean  Green’s  function  is  the  sum  of  two  effective  Green’s  functions.  Finally,  modal 
attentuations  for  the  mean  field  will  be  discussed  here.  It  will  be  shown  that  in  the  case  of 
Dirichlet  boundary  conditions,  in  the  small  roughness  limit,  the  present  method  reproduces  the 
perturbation  results  of  Bass  and  F\iks  [2].  Bass  and  Fuks  also  show  that  when  a  waveguide 
is  not  wide  compared  to  an  acoustic  wavelength,  then  modal  attenuations  arising  from  surface 
roughness  using  half-  space  scattering  amplitudes,  as  in  the  work  of  Kuperman  and  Ingenito 
[3],  are  considerably  different  that  those  calculated  using  renormalized  reflection  coefficents.  It 
should  be  noted  that  Voronovich  has  treated  the  problem  of  waveguides  with  a  single  rough 
boundary  [4)  in  a  somewhat  different  manner.  The  results  presented  here  for  effective  reflection 
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coefficients  are  implicit  in  his  work. 

The  discussion  here  concerns  only  constant  sound  speed  profiles  for  which  up-  and  down¬ 
going  plane  waves  can  be  used.  In  section  II,  scattering  states  will  be  introduced  and  will  be  used 
to  construct  the  held  of  a  point  source  for  each  realization  of  an  ensemble  of  rough  interfaces. 
In  section  III  stochastic  equations  for  the  Green’s  function  will  be  discussed.  In  particular,  a 
self-energy  operator  will  be  defined.  In  section  IV  the  mean  Green’s  function  for  the  case  of 
statistically  homogeneous  roughness  will  be  developed.  Section  V  discusses  modal  attenuations. 
A  summary  of  this  work,  and  an  outline  of  future  work  will  be  found  in  section  VI. 

II  Scattering  states  and  the  construction  of  the  Green’s  func¬ 
tion 

Consider  a  homogeneous  fluid  bounded  by  two  rough  interfaces  as  in  Fig.  1.  These  interfaces 
are  parallel  in  the  mean  and  may  or  may  not  be  penetrable.  The  formalism  to  be  presented  here 
is  not  tied  to  any  particular  boundary  conditions.  However,  assume  that  half-space  scattering 
amplitudes  for  each  of  these  surfaces  considered  separately  are  known:  r'^(K|Q)  for  an  upgoing 
plane  wave  whose  wave  vector  has  a  horizontal  projection  Q  scattering  at  the  upper  interface 
into  a  downgoing  plane  wave  whose  wave  vector  has  horizontal  component  K,  and  r~(K|Q)  for 
a  downgoing  plane  wave  described  by  the  horizontal  wave  vector  Q  scattering  into  an  upgoing 
wave  with  horizontal  wave  vector  K.  In  practice,  one  will  have  only  approximations  for  T*. 
How  can  one  construct  the  Green’s  function  of  the  acoustic  field  using  these  amplitudes?  The 
idea  here  is  that  a  lot  of  physics  has  gone  into  developing  approximations  of  the  scattering 
amplitudes.  Can  this  physics  be  incorporated  directly  into  the  two  interface  problem? 

The  problem  here  is  to  find  the  Green’s  function  which  satisfies 

(V?  -|-u;Vc^)G(r,r,)  =  -«{R  -  R,)«(z  -  z,)  (1) 


where 

r  =  (R,z),ro  =  (Ro.^o). 


In  general,  lower  case  vectors  will  denote  vectors  in  3  dimensions  with  their  horizontal  projections 
given  by  the  corresponding  upper  case  letter.  As  usual,  u  is  the  acoustic  frequency  and  c  is  the 
speed  of  sound  between  the  interfaces.  It  will  prove  convenient  to  work  with  the  double  sided 
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Fourier  transform  of  G, 

6(K,z;Q,Zo)s  Jex^(-iK‘K  +  iq‘Ra)G{K,z]Ro,Zo)dB^-  (2) 

(In  Ref.[l]  only  one-sided  transforms  were  used.  Two-side  transforms  make  the  equations  below 
more  symmetric,  but  they  are  not  absolutely  necesswy.)  The  two-sided  transform  satisfies  the 
separated  wave  equation 

a^G(K,  z;  q,Zo)  +  (wVc*  -  A:*)G(K,  z;  Q,  z„)  =  -{2r)Hiz  -  Zo)SiK  -  Q).  (3) 

At  the  rough  bounding  surfaces  5^,  G  must  respect  whatever  boundary  conditions  are  imposed 
there.  As  a  result,  G  will  not  be  diagonal  in  horizontal  wave  vector.  In  order  to  respect 
these  boundary  conditions,  whatever  they  are,  G  and  G  will  be  constructed  from  two  sets  of 
scattering  states  XqCr)  which  obey  boundary  conditions  on  5^  respectively.  Away  from  the 
boundaries,  these  states  can  be  written  as  superpostions  of  up-  and  down-going  plane  waves 
using  the  associated  half-space  scattering  amplitudes.  By  staying  away  from  the  boundaries, 
the  limits  of  the  Rayleigh  hypothesis  will  not  be  encountered.  In  contrast  to  Ref.[l],  here  the 
Fourier  transforms  of  yiU  be  used: 

X*(K,  Q,  z)  =  exp(±ik,z)«(K  -  Q)  exp(q:i*,z)r*(K|Q).  (4) 

Here  k,  is  the  vertical  component  of  a  wavevector  with  horizontal  projection  K 


It  was  shown  in  Ref.[5]  that  the  scattering  states  are  sufficient  to  compute  half-space  Green’s 
functions.  It  will  be  assumed  that  they  are  also  sufficient  to  compute  the  Green’s  function 
between  two  interfaces.  In  order  to  respect  the  upper  boundary  conditions,  the  Green’s  function 
(j  will  be  constructed  as  a  continuous  superpostion  of  the  states  x'*'  when  z  >  z^  and  as  a  different 
superpostion  of  the  scattering  states  x"  when  z  <  Zo.  The  coefficients  of  the  superpositions  must 
depend  on  Zo  and  the  ‘inddent”  wave  vector  Q.  Hence  G  will  be  written 


6(K,r;Q,Zo) 


/  x+(K,  P,  z)a+(P,  Q,  Zo)dP  If  z>z, 

Jx'iK  p,  ^^)o*(P,  Qi  ^o)dP  if  z  <  Zo 


(5) 


The  a’s  can  be  determined  from  jump  conditions  across  the  plane  contadning  the  source, 
z  =  Zo.  must  be  continuous  across  this  plane  and  the  vertical  derivative  of  G  must  have 
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a  jump  discontmuity  to  produce  the  right  hand  side  of  Eq.  3.  More  explicitly  these  jump 
conditions  are 

G(K,z+;Q,z«)-a(K,z;;Q,Zo)  =  0  (6) 

11  ( 

^d.G(K,z;Q,Zo)|„,+  -^a.G(K,^;Q,Zo)U-  =  (7) 

These  conditions  can  be  written  as  a  matrix-integral  equation  for  the  unknown  coefficients  a. 
Define  elements  of  a  matrix  of  integral  kernels,  M  by 

Mu  (K,P,r„)  =  x'"(K,P,z<,)  (8) 

Mt.2(K,P,^„)  =  -x-(K,P,Zo)  (9) 

Mu(K,P,Zo)  =  (l/ik,)d,x*(K,P,Zo)  (10) 

Mu(K,P,z„)  =  (-l/i*,)d,x-(K,P,z,).  (11) 

The  jump  conditions  across  Zo  can  now  be  written  as 

/  *^u(K.P,^)  *f..,(K,P,.,)  \  /  onp.q,.:)  \  J  «  V  (.2) 

J  ^  M,.,(K,P,j.)  M,,a(K,P,r,)  )  \  o-(P,Q,J,)  )  -‘?g-«(K-Q)  ) 

This  equation  is  essentially  the  Fourier  transform  with  respect  to  the  horizontal  position  of  the 
source,  Ro  of  Eq.  14  in  Ref.[l].  To  condense  the  notation  write  this  equation  as 


with  formal  solution 


which  can  be  written  out  as 


M(zo)a(z„)  =  s, 
q(zo)  =  M~^(zo)s. 


(13) 

(14) 


/  onp,Q.zo)\  r^(  Mi:j(p,K, 
\a-(P,Q,Zo)  )'  J 


zo)  Af,-](P,K,z,)  \  /  0 

z,)  Mll{P,K,z,)  )  I  -Igi«(K-Q) 


(15) 

Note  that  Mu(P,K,Zo)  does  not  denote  the  inverse  of  Mi, i-  Rather  is  the(l,  1)  (operator) 

element  of  M“^. 

The  coefficients  a*  are  now  readily  expressed  in  terms  of  M^  ™  terms  of  M^'j- 

The  resulting  expression  for  the  Green’s  function  G  is 

,)dPl-(2ff)V(*9*)]  if  z  > 
,Zo)dP[-(2x)V(tq,)l  if  z  <  Zo 


C(K,a; <,...)=  I 
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To  this  point  everything  presented  has  been  equivalent,  but  perhaps  more  explicit,  to  the  devel¬ 
opment  in  Ref.[l].  For  flat  surfaces,  the  operator  M  is  diagonal  in  wavenumber  so  that  inversion 
is  simply  the  inversion  of  a  2  x  2  ordinary  matrix.  As  shown  in  Ref.[ll,  the  resulting  Green’s 
function  is  the  same  as  that  given  by  Brekhovskikh  for  flat  interfaces. 


Ill  Finding  the  Mean  Green’s  Function 

The  Green's  function  is  random  because  the  scattering  amplitudes  are  random.  The  primary 
difficulty  in  determining,  for  example,  the  moments  of  G  is  in  finding  the  moments  of  M~*.  In 
Ref.[l]  this  was  done  by  the  method  of  smoothing  which  gives  results  which  only  involve  second 
moments  of  the  fluctuation.  Here  the  problem  of  treating  is  a^idressed  with  methods  of 
quantum  field  theory.  The  idea  is  that  if  one  performs  a  perturbation  expansion  of  am  operator 
about  a  free-space  operator,  long-range  interau:tions  cam  cause  difficulties  in  the  covergence  of 
the  expansion,  amd  even  place  the  existence  of  some  of  the  terms  in  the  expamsion  in  doubt. 
However,  intuitively,  one  supposes  that  prop^ation  between  scattering  events  is  damped  at  long 
ranges  because  a  third  (higher  order)  event  may  come  between  the  two  events  in  question.  One 
might  suppose  that  propagation  between  scattering  events  ought  to  be  described  using  a  mean 
propaigator,  which  is  damped  because  of  the  randomness  of  the  medium.  Of  course  the  mean 
propagator  is  what  one  is  trying  to  And,  so  the  resulting  formalism  will  present  a  problem  in 
self-consistency.  The  formal  aspects  of  such  a  proceedure  will  now  be  presented.  The  operator 
algebra  is  straight  forward,  even  if  the  motivation  may  be  somewhat  obscure.  A  more  thorough 
presentation  of  this  algebra  is  presented  in  the  appendix  of  Ref.[6].  In  contrast  to  the  smoothing 
formalism,  the  algebra  presented  below  results  in  formal  expressions  which  include  all  orders  of 
the  fluctuations. 

Denote  the  inverse  of  Af  by 


g{Zo)  s  M-^{Zo),  OijiZo)  = 


(17) 


The  explanation  of  this  notation  is  that  M~^  is  the  inverse  of  a  linear  operator,  as  is  a  Green’s 
function.  In  order  not  to  confuse  this  inverse  with  the  Green’s  function  G,  the  inverse  of  the 
Helmholtz  operator,  a  calligraphic  Q  is  used. 

Let 


M  =<M  >  +^M. 


(18) 
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The  angular  brackets  denote  an  ensemble  average  over  the  realizations  of  scattering  amplitudes 
Fluctuations  about  the  mean  are  denoted  by  AM.  Just  as  the  inverse  of  M  is  written 

1 


Q  = 


so  the  mean  of  G  is  written 


<(?>  = 


<M>+AM’ 
1 


(19) 


(20) 


<  M  >  +<  E  >■ 

<  E  >is  called  the  average  self-energy  (or  mass)  operator  or  simply  the  self-energy.  The  ter¬ 
minology  comes  from  quantum  field  theory.  Use  Go  to  denote  the  inverse  of  the  average  of 
M: 

I 


ao  = 


(21) 


"  <M>' 

Following  Ref.[6]  it  is  simple  to  show  from  Eqs.  20  and  21  that  the  mean  of  M~^  and  the  inverse 


of  <  M  >  are  related  by 

<G  >  =Go-Qo<'Z><G>> 

This  is  Dyson’s  equation,  eq.  A14  of  Ref.  6.  Fluctuations  of  G  are  found  from 


1 
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^  ^^^~G7^  +  AM  G7' +  <'!:>' 

which  using  simple  operator  algebra  can  be  written  as 

G-<G>  =  -a(AM  -  <  E  >)<  a  >  =  -<  P  >(AM  -  <  E  >)P. 

Finally,  G  can  be  expressed  entirely  in  terms  of  <  p  >  and  a  scattering  operator  T : 

G  =  <G>  +  <G>T<G>, 

where  the  scattering  operator  is  ^ven  by 

r=  -(AM-<  E  >)- 


1 


(22) 


(23) 


(24) 


(25) 


(26) 


H-<P>(AM-  <  E  >)■ 

Expansion  of  the  operator  T  shows  that  the  “propagation”  between  “scattering  events,”  (AM  - 
<  E  >),  occurs  via  the  mean  of  the  inverse  of  M,  <  p  >,  not  via  the  inverse  of  the  mean  of  M, 
Go. 

Since  <p  —  <p>>=:0,  it  follows  that  the  mean  of  the  scattering  operator  T,  vanishes. 


<  r  >=  0. 


(27) 


6 


The  vanishing  of  the  mean  of  the  scattering  operator  T  can  be  used  to  develop  an  expansion  of 
the  self  energy  in  powers  of  the  fluctnations  AM.  Write 

<  E  >  =  +  <^2  +  . . .  (28) 

where  the  <ri  is  the  contribution  to  <  E  >  from  terms  tth  order  in  powers  of  AM.  Equation  27 
and  the  definition  of  the  scattering  operator  T  in  Eq.  26  imply  that 


<To  =  0  (29) 

<ri  =  0  (30) 

ffs  =  -  <  AM<  g  >AM  >  (31) 

ers  =  <AM<g>AM<G>AM>  (32) 

<r4  =  -  <  AM<  g  >AM<  g  >AM<  g  >AM  > 

+<r2<  g  >0-1-  <  AM<  g  ><t2<  g  >am  >  (33) 


For  Gaussian  fluctuations,  ^3  will  vanish.  The  subtractions  in  show  that  it  is  represented  by 
a  skeleton  diagram,  or  an  irreducible  diagram  [6]. 

Substituting  the  expression  for  g  in  Eq.25  in  the  identities  <  MM~^  >=<  M~^M  >=  1 
shows  that  <  G  >  can  be  written  as 

^  ^  ^  ^  <M>  +  <  AM<  g>r>  ^  <  M  >  ^  <T<g  >AM  > 

Comparing  with  the  definition  of  the  self  energy  shows 

<  E  >=<  AM<  g  >T  >=<  T<G  >AM  >  .  (35) 


In  this  equation  <g  >  on  the  right  hand  side  depends  on  <  E  >,  and  T  can’t  be  found  until 
<  ^  >  is  known.  Hence,  although  this  result  is  exact,  it  is  purely  formal.  To  make  progress, 
the  quantities  on  the  right  must  be  approximated.  It  is  to  be  hoped  that  approximation  of  the 
operators  on  the  right  will  be  a  more  effective  proceedure  than  approximating  <g  >  directly. 

It  is  to  be  understood  here  that  all  operators  depend  on  Zo.  In  fact  the  fluctuation  AM  has 
a  simple  dependence  on  Zo. 


AAf(K,Q,Zo)  =  exp(-tfc,Zo)Ar+(K|Q) 


( 


I 

-1 


0 

0 


»exp(+iik.Zo)AT-(K|Q) 


0  1 
0  1 


(36) 
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Thus,  with  summation  over  repeated  indices  and  integration  over  repeated  wavevectors  under¬ 
stood,  the  self-energy  can  be  written  as 


<E(P,Q,^o)  >= 


where 


exp(-ip,z<,) 


( 


A(P,Q, 


A(P,Q,Zo)i.2 
-A(P,  Q,Zo)i,2 


-  exp(-f  tp,Zo) 


( 


A(P  ,Q,  •20)2,1 
A(P ,  Q,  20)2^1 


A(P,Q,  20)2,2 
A(P,Q,  20)2,2 


(37) 


A(P.  Q.  2o)i.i  =  <  Ar+(P|K)<  g  >tAK,  K',  2o)r, •i(K',  Q,  Zo)  >  (38) 

A(P,  Q,  :2o)i,2  =  <  AT-*-(PlK)<  g  K',  Zo)TjAK\  Q,  2o)  >  (39) 

A(P,Q,  20)2,1  =  <  Ar-(P|K)<  g  >2j(K,K',Zo)T>.i(K',Q,Zo)  >  (40) 

A(P,  Q,  zo)2,2  =  <  Ar-(P|K)<  g  >2j(k, k', 2o)r,. 2(K',  q,  zo)  > .  (4i) 

Even  though  the  averages  here  appear  to  depend  on  the  source  level,  Zo,  they,  in  fact,  are 
independent  of  Zo,  as  is  shown  in  the  Appendix.  The  functions  A  will  turn  out  to  be  those  parts 
of  the  effective  reflection  coefficients  not  contained  in  the  mean  half-space  scattering  amplitudes. 

These  formal  expressions  for  g  allow  computation  of  the  averaged  Green’s  function  for  the 
Helmholtz  equation  given  in  Eq.l6.  For  example,  for  the  case  z  >  Zo 

<  G(K,z;Q,z<,)  >=  [-(2)r)V(i9,)l  j  <  Mi.i(K,P,^)  X  ai,2(P,Q,2„)  > 

-I-  <  Ailfi.i(K,  P,  zr)Aai.2(P,  Q,  2,)  >  dP.  (42) 

The  first  product  of  averages  is  easily  expressed  in  terms  of  <  E  >,  and  therefore  A’s,  and  the 
averaged  half-space  scattering  amplitudes,  <T^  >.  The  average  of  the  product  of  fluctuations 
in  Eq.  42  appears  more  problematic,  but  it  too  can  be  expressed  in  terms  of  the  A's,  as  follows. 
The  fluctuation  of  Mi,i  is 

AAfi,i(K,  P,  z)  =  exp(-*it,z)Ar+(K|P).  (43) 

The  fluctuation  of  g  from  eq.25  is  just  <  g  >T<  g  >.  The  average  of  the  product  of  these 
fluctuations  is  (again  with  summations  and  integrations  understood  for  repeated  indices  and 
wavevectors) 
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<  AAfi.i(K,P,z)Aai.i(P,Q,*o)  >= 

exp(-ik,z)  <  Ar+(K|P)  <  ff(P,P',z,)ij  >  T(P',P",Zoh,i  ><  ^0)1.2  > 

=  exp( -ikjz) A(K,P)ij  <  C?(P,Q,Zo)j,j  >  .  (44) 

This  general  expression  for  the  average  of  the  products  of  the  fluctuations  was  missed  in  ref.[l]- 
The  point  here  is  that  the  mean  Green’s  function  for  the  Helmholtz  equation,  G,  is  now  entirely 
expressible  in  terms  of  the  averaged  half-space  scattering  aonplitudes  <  >  and  in  terms  of  the 

functions  A,j.  When  there  is  statistical  homogeneity  it  will  be  shown  that  G  can  be  expressed 
as  the  Green’s  function  for  a  duct  with  flat  interfaces  and  effective  reflection  coefficients  which 
are  combinations  of  <  >  and  the  A’s.  In  terms  of  the  A’s  and  the  mean  of  G  the  Green's 

function  for  the  Helmholtz  equation  is  now  given  by 

{G(K,z;Q,z„)>  = 

/  [exp(+ifc.^)«(K  -P)  +  exp(-tib,z){<T+(K|P)>  +  Ai.,(K,P)}]  (j;(P,Q,^<,)i.2) 

•9*  J 

+  exp(-ifc,z)Ai.a(K,P)(a(P,Q,Zo)2.2)dP,  (45) 

for  the  case  z  >  Zg.  For  the  case  z  <  Zg,  {G)  is  given  by 

(G(K,z;Q,z„))  = 

/  [exp(-i*,z)rf(K  -  P)  -^exp(-^i*,z){(r-(K|P))  -h  A2.2(K,P)}]  <a(P,Q,  ^„)2.2) 

+  exp(-Hi*,^)A2.,(K,P){a(P,Q,«.)i.2)dP.  (46) 
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rV  Statistically  Homogeneous  Roughness 


In  this  this  section  the  results  of  the  previous  section  are  specialized  to  the  case  of  statisticaUy 
homogeneous  surface  roughness.  This  case  is  tractable  because  when  the  statistics  are  homoge¬ 
neous,  the  inverse  of  the  mean  Green’s  function  <  Q  >“*  =  <  Af  >  -I-  <  E  >  becomes  diagonal 
in  wavenumber.  Finding  <  Q  >  itself  then  requires  only  finding  the  inverse  of  a  2  x  2  matrix. 

As  discussed  in  the  appendix  to  Ref.{l],  statistically  homogeneous  roughness  means  that  if 
z  =  h(R)  is  a  member  of  a  ensemble  of  surfaces,  then  z  =  Z>aA(R)  =  A(R  -  a)  is  an  equally 
probable  member  of  that  ensemble,  where  a  is  any  constant  translation  vector  in  the  x-y  plane. 
In  Ref.[l]  it  is  also  argued  that  scattering  amplitudes  are  functionals  of  the  surface  function  h(R). 
(Here  there  are  two  surface  functions,  one  for  the  upper  surface  and  one  for  the  lower  surface. 
When  necessary  h  will  be  understood  to  include  both  surfaces,  and  translations  will  apply  to 
both  simultaneously.)  If  this  functional  dependence  is  made  explicit,  as  in  T(K|Q,[h]),  then 
simple  physical  arguments  show  that  the  scattering  amplitudes  for  translated  and  untranslated 
surfaces,  Dgih  and  h  are  related  by 

T*(K|Q,  (/»])  =  exp{*(K  -  Q)  •  a)T*(K|Q,  [D^h]).  (47) 

Applying  the  same  translation  to  both  the  upper  and  lower  surfaces  shows  that  the  matrix 
operator  M  behaves  similarly: 

M(K,  Q,  [AD  =  exp(t(K  -  Q)  •  a)Af(K,  Q,  [D^h]).  (48) 

Since  Q  is  the  inverse  of  Af,  it  too  behaves  this  way  under  horizontal  translations  of  the 
waveguide: 

a(K,  Q,  [h])  =  exp(i(K  -  Q)  •  a)^(K,  Q,  [Dt,h]).  (49) 

Because  of  the  assumed  statistical  homogeneity,  it  follows  that 

<  a(K,  Q,  Zo,  [h])  >=<  a(K, Q, z„,  [£>a*])  >,  (50) 

for  aU  translations  a  and  all  wavevectors  K  and  Q.  This  means  that  <  ^(K,Q,  Zo.[A])  >  must 
be  proportional  to  a  delta  function  in  K  -  Q,  i.e.  it  must  be  diagonal  in  wavevectors,  and  that 
it  can  be  written  as 

<  a(K,Q,z„,[A])  >=  6{K  -  Q)ff(K,Zo).  (51) 
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When  eq.  25  is  used  to  compare  ^’s  for  translated  and  untranslated  surfaces,  it  follows  from 
eqs.  49  and  51  that  T  has  the  same  properties  under  translation: 

r(K,Q,ro,(/i])  =  exp(t(K  -  Q)  •  a)T(K,Q,  (52) 

The  matrix  elements  of  the  average  self  energy  S,  the  A’s  (treated  generically)  satisfy 

A  =<  AT(K|P,[/i])i/(P,r<,)T(P,Q,ro,[/»l  >  = 
exp[*(K  -  Q)  •  a]  <  AT(K|P,  Da[*))p(P.  «o)r(P,  Q,  r,,,  >  •  (53) 

The  averages  are  equal  for  all  translations  and  it  follows  again  that  the  A’s  are  diagonal  in 
wavevector: 

A.j(K,  Q)  =  S(K  -  Q)Aij(K).  (54) 

Obviously  the  averaged  self-energy  <  £  >  is  also  diagonal  in  wavevector  and  can  be  expressed  in 
terms  of  the  A’s.  The  same  arguments  show  that  the  mean  scattering  ampltudes  can  be  written 

<T*(K|Q)>=^(K-Q)t*(K).  (55) 

It  is  now  possible  to  express  the  mean  Green’s  function  g  in  terms  of  the  A’s.  The  function 
g  is  found  from  the  inverse  of 


<  M (K,  Q,  ^o)  >  +  <  S(K,  Q,  z„)  >=  «(K  -  Q)  X 

(-)-exp(-iA:,z<, )(«■*■  +  Ai.i)  -I-  exp(tA:,z„)(l  -  Aj.i),  -exp(-»*^Zo)(l  -  A1.2)  -  exp(-|-tA:*z<,)(r  +  A2.2) 

-exp(-iA:jZo)(t+  +  Ai,i)  exp(tA,2o)(i  -  A2,i),  -f-exp(-iAtZ<,)(l  -  Ai,2)  -  exp(-|-i*:2Zo)(<~  +  A2.2) 

(56) 

and  is  given  by 

g{K,Zo)  =  2D{K) 

(-)-exp(-tfc,z<,)(l  -  Ai,2)  -  exp(tfcz2o)(^~  +  ‘^2.2).  +exp(-tfczZo)(l  -  Ai,2)  +  exp(ik^Zo)(t~  +  A2,2)  ^ 

-  exp(tA:,Zo)(l  -  Aa,i)  +  exp(-ifcz2o)(t'''  +  exp(ifc,2o)(l  -  ■^2,1)  +  exp{-iktZo){f^  +  Ai,i)+  ^ 

(57) 

It  is  understood  here  that  t*  =  t*(K)  and  =  Xij(K).  The  determinant  2D  is  given  by 


2i?(K)  =  2[(1  -  Ai,2(K))(l  -  A2,i(K))  -  (t+(K)  -I-  Ai,i(K))(t-(K)  -I-  A2.2(K))1  (58) 
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These  results  can  now  be  used  to  write  the  mean  Green's  function  of  the  Helmholtz  equation, 
G,  in  terms  of  Xij  and  Combining  Eqs.  42,  44,  51  amd  54  gives 

<  G(K,  z,  Q,  z,)  >=  -  Q)G(K,  z,  z,).  (59) 

For  z  >  Zo  the  reduced  Green’s  function  is 

G(K,z,Zo)  =  (exp(i*:,z)  +  exp(-*fc,2)(t+(K)  +  Ai,i(K))]si,2(K,  Zo)+ 

exp(-iA:*z)Ai.2(K)p2,2(K,  Zo),  (60) 

while  for  z  <  Zo  it  is  given  by 

G(K,z,Zo)  =  [exp(-iktz)  -h  exp(iktz}(r(K)  +  A2,2(K))]p2,2(K,  Zo)+ 

exp(  +  ifc,z)A2,i(K)pi,2(K,Zo).  (61) 

Equation  57  gives  expressions  for  the  gij  that  appear  here  in  terms  of  t's  and  A’s. 

Reciprocity  in  the  form 


G’(R>,  z;  Rq,  z^)  —  ^(1101  ^0?  ^)» 

(62) 

implies 

G(K,z;Q,z,)  =  G(-Q,Zo;-K,z), 

(63) 

and 

G(K,  z,  Zo)  =  G(-K,  Zo,  z). 

(64) 

Comparing  coefficients  of  exp(ik,(z  -  Zo))  in  these  expressions  shows  that 

(1  -  A2.i(K))/I>(K)  =  (1  -  Ai,2(-K))/Z?(-K). 

(65) 

If  isotropy  is  assumed  in  the  form 

A.  *'.)  =  M-K), 

and 

t*(K)  =  <*(-K) 

then  Eq.  65  implies  that 

Ai.2(K)  =  A2.i{K). 

Then,  assuming  isotropy,  it  is  possible  to  combine  the  parts  of  G(K,  z,  Zo)  to  show  what  is  the 
principal  result  of  this  work, 
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-  V  (<•*■■> +.-*->V+,(K)|(.-“-'<+.+'*-<n7/(K)|, 

'  2(1  -  A,j(K))(1  -  n//(K)V,7/(K)) 

Ai.2(K)exp(-iA;z|-g  -  2o\)  ,  . 

2(1-Ax.2(K))  ^  ^ 

In  this  equation,  z>  is  the  greater  of  z  and  z<,  and  z<  is  the  lesser  of  z  and  z,,.  The  effective 
reflection  coefficients,  are  given  by 

■  (1-Ax,2) 

Kn  =  (68) 

When  there  is  no  correlation  between  the  fluctuations  in  the  scattering  amplitudes  at  the  upper 
and  lower  interfaces,  Ai,2  will  vanish,  at  least  to  second  order  in  the  fluctuations.  Then  G  is 
exactly  of  the  form  described  by  Brekhovskikh  [7]  for  flat  surfaces  with  reflection  coefficients 
This  is  also  the  form  found  in  ref.[l]  using  the  method  of  smoothing.  Here  the  existence  of 
effective  reflection  coefficients  is  demonstrated  to  all  orders  in  the  fluctuations  of  the  scattering 
amplitudes. 

The  mean  field  will  be  d'escribed  by  normal  modes  which  are  determined  from  the  poles 
of  G(K).  Apparently  these  poles  can  arise  from  the  zeros  of  1  -  Ai,2  or  from  the  zeros  of 
1  -  V^j{K)V~jf{K).  However,  as  Ai,2  -►  1,  -►  oo.  Keeping  track  of  dominant  terms  shows 

that  in  fact,  G  is  not  singular  as  Ai,2  — »  1.  Thus  the  normal  modes  are  determined  only  by  the 
zeros  of  1  -  These  zeros  occur  for  complex  K  and  the  imaginary  parts  of  K 

determine  modal  attenuations. 

In  the  next  section  the  A’s  will  be  discussed  in  more  detail  in  the  lowest  order  approximation, 
when  only  second  order  terms  in  the  fluctuations  of  scattering  amplitudes  are  kept.  It  turns  out 
that  the  effective  reflection  coefficients  describe  all  processes  in  which  a  wavevector  is  forward 
scattered.  The  mean  half-space  scattering  amplitude  includes  only  processes  than  involve  a 
single  interface.  When  a  second  interface  is  present,  there  are  additional  processes  by  which  a 
wave  may  be  forward  scattered.  These  processes  are  contained  in  the  Ai,i  and  A2,2  additions 
to  the  mean  half-space  scattering  amplitudes,  to  form  the  effective  reflection  coefficients  for  the 
mean  field  between  two  interfaces. 


V  Second  order  calculations 


In  this  section  the  simplest  approximation  for  the  effective  reflection  coefficients  will  be  de¬ 
veloped.  In  effect,  this  means  finding  the  average  self-energy.  It  wUl  be  shown  that  if  this 
result  is  further  approximated  by  using  pertubative  approximations  for  the  half-space  scattering 
amplitudes,  then  the  approximation  of  Bass  and  Fuks  [2]  is  obtained. 

To  find  the  self-energy,  an  approximation  for  the  scattering  operator  T  is  required.  Of  course, 
to  find  T,  <  £  >  is  required,  but  one  has  to  start  somewhere.  Equations  28-33  show  that  the 
self-energy  will  be  at  least  second  order  in  fluctuations  of  the  half-space  scattering  amplitudes. 
This  means  that  a  first  order  approximation  for  T  (see  £q.26)  is  simply 

r  «  ~AA/,  (69) 

and  this  is  the  approximation  that  will  be  used  throughout  this  section.  This  approximation  does 
not  require  that  AT  be  expressed  in  powers  of  the  surface  roughness.  One  could  approximate 
AT,  for  example,  by  the  first-order  small  slope  approximation. 

In  any  event,  if  for  the  purposes  of  this  section  we  assume  Eq.  69,  then 

Ai,i(P,Q)  =  -J dKffiAKz,)  <  AT^(PIK)AM,j(K,Q,Zo)  >,  (70) 

where  AT"*"  is  to  be  used  when  i  =  1  and  AT~  when  i  =  2.  For  example,  from  Eq.36  it  follows 
that 

Ai.i(P,Q)  =  j dK[-pi.i(K,Zo)  +  9i.j(K,Zo)lexp(-tM„)<  AT+(P|K)AT+(K|Q)>.  (71) 

From  eq.  57,  the  sum  of  the  p’s  here  is 

_  /XT  _  \  _  /v  _ /-I.  _  ^(<■(K)  + A2,2(K)) 

~  Zo)  -f- 

SO  that  Ai,i(P,  Q)  is  given  by 

Ai.i(P,  Q)  =  y  <  AT+(P1K)AT+(K1Q)  >  (73) 

Note  how  the  exponential  dependence  on  the  source  level  has  dissappeared.  Statistical  homo¬ 
geneity  implies  that  the  average  of  the  fluctuations  of  the  scattering  amplitudes  can  be  written 
as 

<  AT=^(P|K)AT*(K|Q)  >=  «(P  -  Q)c=^'=*=(P,  K),  (74) 
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so  that 


Au(P,Q)  =  «(P-Q)  I 

=  ^(P-Q)Au(P).  (75) 

One  can  show  similarly  that 

Aj.2(P)  =  j  (76) 

Ai.j(P)  =  J  dKi^—^^:^c+<-{P,K).  (77) 

These  three  equations  form  a  coupled  set  of  integral  equations;  one  can’t  find  Ai,i  without 
knowing  Aj^ji  fof  example.  If  there  is  no  correlation  between  the  fluctuations  in  scattering 
amplitudes  at  the  upper  and  lower  interfaces,  then  e'^'~  will  vanish  and  Ai,2  along  with  it. 
Furthermore,  Ai,i  and  Aj,}  are  second  order  in  the  fluctuations  and  should  be  small  compared 
to  t^.  Thus,  the  A’s  are  further  approximated  by 

A,,,(P)  =  »  K)  (78) 

*J.J(P)  =  J  (79) 

To  this  same  order  the  determinant  i)(K)  should  be  written  as 

Z?(K)  =  1  -  «+(K)r  (K).  (80) 

There  might  be  zeros  in  I1(K)  near  the  effective  modal  wave  numbers.  However,  since  are 
the  mean  half-space  scattering  amplitudes,  the  zeros  of  J?(K)  should  occur  for  complex  AT. 

One  can  expect  that  Eqs.  78  and  79  provide  reasonable  approximations  if  the  A’s  determined 
by  these  equations  are  small  compared  to  f**  and  t~.  It  is  also  necessary  that  the  zeros  of  D(K) 
not  be  shifted  significantly  by  the  use  of  the  approximation  in  Eq.  80.  Without  doing  detailed 
calculations  it  is  difficult  to  determine  when  the  A’s  will  be  small.  However,  consider  the  following 
“toy”  problem.  Suppose  that  the  fluctuations  in  scattering  amplitudes  are  sharply  peaked  in 
the  forward  direction,  corresponding  to  an  ensemble  of  nearly  flat  surfaces.  In  fact  suppose  that 

c+'+(P,K)  =  <t+«{P-K)  (81) 

c-’-(P,K)  =  <r-tf(P-K).  (82) 
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Then  Eqs.  78  and  79  |dve 


At.i(P)  = 
Aj.2(P)  = 


t-(P)a+ 

D(P) 

t'»-(P)<T- 

D(P) 


The  approximations  in  Eqs.  78  and  79  will  be  seif-consistent  provided 

r(P)<r+ 


Z?(P) 

D(P) 

These  two  inequalities  in  turn  reqture  that 


<  f"(P) 

<  r(P). 


(83) 

(84) 


(85) 

(86) 


Z)(P)>  ^ 


(87) 


If  .the  zeros  of  D  occur  for  real  values  of  P  this  inequality  will  be  violated.  Otherwise,  it  would 


seem  possible  to  satisfy  this  construnt  if  the  strengths  of  the  fluctuations  are  sufficiently  small. 


The  corresponding  approximations  for  the  effective  reflection  coefficients  now  become 

V3,(P)  =  +  (88) 

K},m  =  ,+(P)  +  |ac^_Jtm_,-.-,p,K).  (89) 


If  there  is  a  small  correlation  between  AT*  and  AT*,  then 


and  the  effective  reflection  coefficients  will  be  reduced  by  1  -  Ai,2.  These  expressions  for  the 
effective  reflection  coeficients  were  derived  in  ref.[l]  by  the  method  of  smoothing.  Here  they  have 
been  shown  to  be  low-order  approximations  for  more  general  expressions.  The  result  for  Ai,2  is 
new  here.  It  is  not  known  yet  how  important  Ai,2  is,  but  it  seems  likely  that  it  is  important  in 
scattering  in  thin  sedimented  layers  where  interfaces  tend  to  follow  one  another. 

The  expressions  for  the  effective  reflection  coefficients  can  be  spedalized  even  further  in  the 
case  of  small  surface  roughness.  This  spedalization  will  be  shown  here  to  compare  the  pres  mt 
results  with  those  of  Bass  and  Fuks  [2]  and  to  illustrate  the  contribution  of  the  self-energy  to 
the  effective  reflection  coefficients. 
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Consider  two  Dirichlet  surfaces  on  which  the  field  vanishes,  and  which  are  described  by 

r  =  if*  +  fi*(R).  (91) 

The  nominal  width  of  the  layer  between  these  sufaces  is  L  =  R+  -  H~  >  0  and  the  mean  values 
of  roughnesses  vanish.  To  second  order  in  the  surface  roughness  the  scattering  amplitude 
for  the  upper  surface  with  a  plane  wave  incident  from  below  is  [8]  (with  exponentials  appearing 
because  the  mean  surface  is  not  at  z  =  0) 

^^(KIP)  as  -  P)  -  -  K) 

+  e'(^*+**)"*  j  dQ2p,?,A+(Q  -  K)fi+(P  -  Q),  (92) 

where  the  Fourier  transform  of  the  surface  roughness  is  defined  by 

(93) 

For  a  plane  wave  incident  from  above  on  the  lower  surface  whose  mean  is  at  R~,  the  scattering 
likewise  is 

r-(K|P)  «  -  P)  +  2ip,e-*<'*+^*)"”fi-(P  -  K) 

+  J  dQ2p,?,A-(Q  -  K)fi'(P  -  Q).  (94) 

The  spectra  of  the  surface  roughnesses,  5*'*,  are  defined  from 

<  /i*(Q  -  K)fi*(P  -  Q)  >=  SiP  -  K)5*’*(P,  Q).  (95) 

Through  second  order  in  the  surface  roughness,  the  mean  half-space  scattering  amplitudes  are 
seen  to  be 

t=*=(K)  =  -I-  2k,  j  dQ?,5*’*(K, Q)}.  (96) 

On  the  other  hand,  the  correlations  of  the  scattering  amplitudes,  are  given  by 

c+’+(K,P)  =  -4*,p,5+’+(K,P)e-*-*’**"^+*’'’'"^  (97) 

c-’-(K,P)  =  -4k,p,5-’-(K,P)e-***‘"‘ (98) 

c+--(K,P)  =  -(-4ib,p,5+*-(K,P)e+’**(^''-"")+‘'’*("''-"“)  (99) 

To  compute  the  self-energies,  Aj^-,  approximate  the  determinant  D  using  the  flat  surface  reflection 
coefficients,  tf 

DiK)  w  1  -  ttt~  -  1  -  =  1  -  e*‘**^.  (100) 
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In  tius  perturbative  approximation,  the  self-energies  are  given  by 


Ai.x(K)  = 

(101) 

Aj,2(K)  = 

(102) 

Ai.2(K)  = 

f  ^2i{kt+PM)L 

(103) 

As  indicated  in  the  previous  section,  the  pole  structure  in  the  integrands  of  these  expressions 
is  capturing  the  fact  that  in  a  layer,  forward  scattering  can  occur  not  only  by  first  scattering 
to  an  intermediate  wave  number  P  and  then  back  to  K  by  a  single  surface,  as  expressed  in 
the  various  contributions  to  <  >,  but  also  that  forward  scattering  can  occur  by  scattering 

into  P,  then  having  any  number  of  unperturbed  forward  scatterings  when  the  energy  bounces 
between  the  interfaces  of  the  layer,  and  then  rescattering  into  the  incident  wavevector  K.  The 
denominators  in  these  expressions  for  A  are  the  phases  acquired  in  bouncing  between  interfaces. 
One  can  see  that  the  self-energy  contribution  to  effective  scattering  amplitudes  will  have  a  much 
stronger  frequency  dependence,  because  of  the  denominators  in  the  integrands,  than  the  half¬ 
space  scattering  amplitudes. .  These  contributions  are  of  the  same  order  of  magnitude  as  the 
perturbations  of  the  flat  surface  mean  scattering  amplitudes. 

In  fact,  t*  and  Ai,i  can  be  added  to  show  how  the  result  of  Bass  and  Fuks  can  be  obtained 
from  this  formalism. 


.  r  1 4. 

t+(K)  -I-  Ai.i(K)  = 


(104) 


This  expression  can  be  written  in  terms  of  the  Green’s  function  between  two  flat  Oirichlet 
surfaces  at  From  Brekhovski  or  from  eq.66  above  this  Green’s  function  is 


Go(K,  z,  Zo) 


(2ir)^  [«’*»*>  - 

-2tfc,  1  - 


(105) 


It  is  now  straight  forward  to  show  that  the  integrand  in  the  effective  reflection  coefficient  can 
be  expressed  in  terms  of  a  second  derivative  of  Gy. 


lim  lim  d,d,«G<,(K,s,Zo)sG"(K) 

*o-»W+  x—H* 

1  1  fZiktL 


(106) 
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The  effective  reflection  coefficient  now  can  be  written  as 

V*„(K)  =  «“■«•{-!  -  j (107) 

Bass  and  Fuks  calculate  an  effective  impedance,  which  is  the  ratio  of  the  mean  field  to  the 
z-derivative  of  the  mean  field.  When  the  effective  reflection  coefficient  is  the  effective 
impedance  m  is 

m(K)  =J^f  rfPp,G"(P)S+*+(K,P).  (108) 

Apart  from  factors  of  k  and  2  which  arise  from  differences  in  normalization  of  the  Green's 
function  and  the  Fourier  transforms,  this  is  precisely  the  impedamce  of  Bass  and  Fuks. 

VI  Summary 

This  work  is  an  extension  of  the  ideas  presented  in  Ref.fl].  The  primary  result  is  Eq.  66 
which  shows  that  the  mean  Green’s  function  in  a  medium  between  two  rough  interfaces  can  be 
expressed  as  a  weighted  sum  of  two  familiar  Green’s  functions.  The  first  is  the  Green’s  function 
in  a  layer  bounded  by  flat  interfaces  with  effective  reflection  coefficients,  The  second  is 
a  free-space  Green’s  function  whose  exponent  has  a  sign  opposite  to  that  of  the  conventional 
free-space  Green’s  function.  The  z-derivative  of  the  combined  Green’s  function  has  the  correct 
jump  discontinuity  across  the  plane  z  =  Zq.  This  result  extends  that  of  Ref.[l]  in  that  it  is 
not  restricted  to  be  second  order  in  fluctuations  of  half-space  scattering  amplitudes.  It  seems 
to  be  a  quite  general  result  that  surface  roughness  can  be  incorporated  into  the  mean  field  by 
effective  reflection  coefficients.  In  addition,  here  the  possibility  of  correlated  surface  roughness  is 
considered.  To  show  that  the  methods  used  here  are  not  unreasonable,  results  were  specialized 
to  first  include  only  second-order  terms  and  second,  to  include  only  terms  second  order  in  surface 
roughness.  In  this  case,  the  results  of  Bass  and  Fuks  were  obtained. 

The  motivation  for  this  inquiry  is  to  find  a  way  to  incorporate  what  is  already  known  about 
half-space  scattering  into  the  waveguide  problem.  One  shouldn’t  have  to  do  the  half-space 
problem  again  in  the  waveguide.  For  example,  with  the  results  presented  here,  one  could  consider 
propagation  in  a  fluid  layer  bounded  by  an  elastic  solid  without  repeating  some  fairly  messy 
computations. 

The  effective  reflection  coefficients  are  not  simply  the  mean  half-space  scattering  amplitudes. 
Instead,  these  effective  reflection  coefficients  account  for  channels  of  forward  scattering  that  are 


aot  avulaUe  when  there  is  only  one  interface.  Modal  wave  numbers  are  found  from  the  dispersion 
rdation 

i-v*,(K)v;7,(K). 

As  discussed  by  Bass  and  Fuks,  the  modal  wave  numbers  can  be  quite  different  from  those  found 
using  only  the  half-space  scattering  amplitudes. 
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Appendix 


It  is  argued  here  that  in  the  self-energy  ^ven  in  Eq.  37,  the  A’s  must  be  independent  of  the 
source  level,  z^.  This  is  important  for  it  means  that  the  effective  reflection  coefficients  for  the 
mean  field  will  be  independent  of  the  source  level,  a  desirable  property.  When  A],3  ^  0  it  appears 
that  this  will  not  be  the  case,  even  though  the  A's  are  independent  of  Zo- 
First  note  that  the  derivative  of  M{z)  satisfies 


dM(z)  _ 
dz  ~ 


1 

0 


Af(z). 


(Al) 


Since  by  definition 


M{z)Miz)-'  =  1, 


(A2) 


the  derivative  of  the  inverse  of  M  satisfies 


dAf-^(z) 

dz 


Averaging  this  equation  shows  that 

d  <  Qjz)  > 
dz 


<Q> 


(A3) 


(A4) 
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Repeating  these  argoments  shows  first  that 


°  '  <c>- 

\l  0  ) 

This  equation  obviously  holds  for  <  Af  >,  and  since  <  Q  >“^=<  M  >  +<  S  >, 

d< 


(A5) 


<  E  >. 


(A6) 


Using  Eq.  37  for  <  E  >  shows  that  A  must  satisfy 


exp(-tp,z)^^^  -  exp(+tp,z)^^^  =  0, 


(A7) 


where  Ai  is  the  first  matrix  in  eq.37  and  A3  the  second. 

Equation  (A6)  can  be  differentiated  again  to  show  that  the  self-energy  satisfies  the  separated 
wave  equation: 

E  >  =  0.  (A8) 

This  means  that  the  A’s  must  be  of  the  form 


Aij  =  Aj -k- Bj  exp{2ip,z)  (A9) 

Asj  =  Cj  +  Djexp{-2iptz)  (AlO) 


with  A,  B,  C,  D  constant.  Then  eq.A7  implies 


Bj  exp(tpxz)  -  Dj  exp(-tp,z)  =  0. 


(All) 


The  only  way  this  can  hold  over  a  continuous  range  of  z  is  for  the  B's  and  D’s  to  vanish,  which 
implies  that  the  A’s  are  independent  of  source  level. 
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Figure  Captions 

Figure  1.  A  sketch  of  the  geometry  used  in  this  paper.  The  z-axis  is  positive  upward.  The  dashed 
line  indicates  an  imaginary  plane  through  the  source.  Between  the  planar  surfaces  shown  here, 
the  fleld  can  be  expanded  in  half-space  plane- wave  scattering  states,  These  states  describe 
scattering  from  a  wavevector  K  into  a  number  of  plane- waves  with  wavevectors  Q. 
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